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Abstract 

We consider the reaction zone that grows between separated regions of dif- 
fusing species A and B that react according to mA + uB — > 0, within the 
framework of the mean- fieldlike reaction-diffusion equations. For distances 
from the centre of the reaction zone much smaller than the diffusion length 
Xd = y/Dt, the particle density profiles are described by the scaling forms 
predicted by a quasistatic approximation, whereas they have a diffusive cut- 
off at a distance of order X^. This cutoff, and the power-law decay of the 
quasistatic profiles, give rise to multiscaling behaviour, with anomalous val- 
ues for the exponents describing the moments of the density and reaction 
profiles. Numerical solutions of the reaction-diffusion equations are in good 
quantitative agreement with the predictions of this theory. 

PACS numbers: 82.40.-g, 82.30.-b, 02.30.Jr, 02.70.Bf 
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I. INTRODUCTION 



The problem of a front that grows between initially separate regions of diffusing species 
A and B that react according to mA + nB has provoked much recent interest . 
Studies have concentrated on the scaling properties of the reaction rate R per unit volume 
and particle density profiles a and b, and on the critical dimension above which the mean- 
fieldlike react ion- diffusion rate equations are valid. In geometries where the only spatial 
variation of the densities is along the x-axis, these equations take the form [|I],Q 

dta = Ddla - mka'^b'', (1) 
dtb = Ddlb - nka"'b'', (2) 

where k is the reaction constant, and here and subsequently the reagents are assumed to 
have the same diffusion constant D. The initial conditions appropriate to this problem are 

a{x, 0) = aQ9{x), (3) 
b{x,0) = boe{~x), (4) 

where Oq and bo are constants, and 9 is the Heavyside function. The quantity u = {a/m—b/n) 
obeys a simple diffusion equation, with solution p 

. = if^-M-if^ + Merff^V (5) 

The reaction is concentrated in the region where the densities of the two species are com- 
parable. The centre x/ of the reaction zone may be defined as the point where u = 0. We 
therefore have 



Xf = 2(Dt)^/2g^f- 



(gp/m - bp/n) 
(oo/m + bo/n) 



(6) 



so if ao/m = bo/n we have Xf = (n.b. in reference 0, x/ was defined as the point of 
maximal reaction, which is not necessarily the same point). 

If one assumes that the penetration of one species into the other is much shallower 
than the diffusion length (Dt)^/^, the reaction between the two species takes place within a 
distance w <^ {DtY^"^ oixf. One expects that, for x <^ (DtY^'^, the profiles will be described 
by the single lengthscale w, leading to the following scaling hypothesis plJ 



a(x,t)=t-^A(^^), (7) 
bix,t)=t-^B(^^), (8) 



t' 

i?(x,t) = r/^0(^j, (9) 

where w ^ t", (p = A'^B"', and (3 = {m + n)7. The number of particles of either species 
arriving at the reaction front is oc which must equal the total reaction rate, so we 
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must have (3 — a = 1/2. Consistency of this scahng ansatz with the equations of motion 
leads to a = (l/2)(m + n- l)/(m + n + 1) H. 

The related case of a front formed by opposing constant diffusion currents J4 = mJ and 
Jb = —nJ of A- and 5-particles imposed at x = —00 and +00 respectively has recently 



been studied [p!0H12|. In this case, the system approaches a steady state where the equations 

iD/m)dla = ka^'b" = {D/n)dlh, (10) 
and boundary conditions may be written in dimensionless form, so that the following scaling 



ansatz is valid for all x |11] 



i?(x) = -0,, - , (11) 
w \wj 

a{x) = JwA,, , (12) 

h{x) = JwBss (-) , (13) 



where w{J, D, k) oc J^*^ and u = {m + n — l)/{m + n + 1). From ([10|), we have d^u = [u = 
{a/m—b/n)], whose solution with these boundary conditions is m = —Jx/D. For (x/w) ^ 1, 
the 5-particles are overwhelmingly in the majority, so one has b = nJx/D+na/m ^ nJx/D, 
and hence A'4(l/) ~ (^ss)™"?/", leading to 

A (v) ~ /l/""^^exp(-(Ty^+("/2)) for m = 1, . 
Ass[y) ~ i .,_(„+2)/(,n-i) for m > 1 ^ ^ 



y 



as y 00, where a is a constant. Similar results hold for B^g by interchanging m and 



n. Within this approach, it is also possible to show ||TT| that the 'mean-field' assumption 
R = ka"^b"' (assumed in all the above equations) is valid for microscopic stochastic systems 
in spatial dimension d > dc = 2/ {m + n — 1). 

For the time- dependent problem, when x ^ w one of the species is overwhelmingly in 
the majority, so |n| ma.x{a/m,b/n), and the profile of the majority particle density is 
~ \x — Xfl/t^^"^ for (Dty^"^ ^ X ^ w. The diffusion current of particles arriving at Xf 
is therefore J ~ t^^/^, and the characteristic timescale on which this current changes is 
{dlog J/dt)~^ oc t. The equilibration time of the front is of order Dw^, so since a < 1/2 
one would expect that the reaction zone has enough time to reach the steady-state profile 
it would have if the current J were constant. One would therefore predict that the results 
of the steady-state problem, and hence the dynamic scaling ansatz, would be applicable to 
the time dependent case for x <C {Dty/"^ |]lT|. 

For m = n = 1 the scaling forms (0^ have been proved rigorously to describe the 
asymptotic behaviour as t — cxo of the reaction-diffusion equations (|I],0) ||^. Experiments 
on real systems, and simulations of microscopic stochastic models, also appear to verify the 
scaling theory and exponents in dimension d > 2 [0-^. For d = 1, there has been some 
controversy as to whether the scaling theory is valid |^,|], but the most recent results 
appear to show that the steady-state results do indeed apply. However, for (m,n) 7^ (1, 1), 
where rigorous mathematical results are not available, the case is much less clear. Numerical 
simulations of microscopic stochastic models in d = 1 are consistent with a scaling ansatz 
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0], but they are of low precision, since reaction events are much rarer than for m = n = 1. 
The fact that at least one of the particle density profiles must decay algebraically (see Eq. 
(p!^)) might invalidate the assumption that reactions take place within a zone of width 

In this paper, we shall give careful arguments to show that the scaling ansatz is indeed 
valid for lengthscales much smaller that ~ t^^"^. We shall then show that the two lengthscales 
in the problem, w and {Dty^'^, together with the power law tails of the steady-state profile, 
give rise to a multiscaling form for the particle density profile, whose moments are described 
by a spectrum of exponents between a and 1/2. We then present numerical solutions of 
the reaction-diffusion equations, and show that the results are in good agreement with the 
theoretical predictions. 



II. VALIDITY OF THE SCALING ANSATZ 

The scaling ansatz can be shown to be exact for the case of a steady-state front formed 
between balancing opposing currents, and so its applicability to the time-dependent case 
relies on the front being formed quasi-statically Naively, one would expect that the 



time for a diffusive system to equilibrate within a region of size ~ t" would be oc t"^", whereas 
the timescale upon which the current J oc t"^/^ changes is (dlog J/dt)~^ oc t, which predicts 
that quasistatic approximation would be valid for lengthscales with v < 1/2. However, 
since some of the density profiles decay algebraically, one might wonder whether the flow 
of particles towards |x| ^ oo necessary to sustain these steady-state profiles might be too 
great for the quasi-static approximation to be valid. 

In this section, we shall show that the quasistatic approximation is internally consistent 
for lengthscales smaller than t^/^, in that it predicts that: (i) the number of particles up to 
a distance ~ t^/^ is always much less than the total particles that have reacted, so that the 
number of particles in the tails is never too much to have a feedback effect on the profiles 
at distances of order t^/^; (ii) the time taken for each part of the particle density tail to 
equilibrate at its quasistatic value is always much less that the characteristic timescale on 
which this value changes. We shall then discuss what this implies about the behaviour of 
the moments of the density and reaction profiles. 

Consider the part of the tail of the A-particle profile a(a;, t) ~ t-T(x/r)-^, where A = 
[n + 2]/[m — 1] (see Eq. (p!4D ), in the region xi < x < X2, with xi oc t"^ and X2 oc 
(a < ei < 62 < 1/2). The current of A-particles at x is 

Ja = -Dd^a ~ (15) 

X 

so the ratio JA(a^2)/-^yi(a^i) = if:~('^2-ei)(i+A) q g^g ^ _^ Almost all of the particles that 
enter at Xi are therefore removed by the reaction, rather than by diffusing out at X2- The 
number of particles in the region is 

{a{xi)xi for A > 1, 

r-^log(f^) forA = l, (16) 
a{x2)x2 for A < 1. 
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The number of particles in the tail can diverge for certain values of A. This could invalidate 
the assumption that the total reaction rate equals the number of particles arriving at the 
origin if this number were found to be larger than the total number of particles oc t^/^ that 
have reacted. However, the total number of particles in the tail up to a lengthscale t^/^, 
found by substituting e2 = 1/2, is found in each of the above cases to be of order less than 
The time taken for A^^ particles to enter the region a;i < x < 0:2 is 

r t^'' for A > 1, 

Na/Ja{Xi) ~ <^ t^ei iQg^ fo^ ^ _ 1^ (17) 

[ i.l+e2-(.2-ei)(l-A) ^ > 1^ 

which is always ^ t since ei < £2 < 1/2. The front therefore has enough time to reach its 
steady-state value for lengthscales smaller than t^/^. 

The scaling ansatz (^^) would therefore appear to be consistent for all lengthscales ~ f^, 
with e < 1/2. For e > 1/2, the density of particles is limited by diffusion, and so we expect 
there to be some kind of exponential cutoff in all of the profiles on such lengthscales. We 
therefore propose the following ansatz for a and b in the limit t ^ 00: 

a{x,t) = ass{x,t)GA ' (18) 

h{x,t) = hss{x,t)GB[^-^y (19) 
where ass = t~'^ Ass{x / 1^) and hgs = t~'^ Bss{x /t^) are the solutions to the steady-state 



equations (|12| , p^ , and GAiy) and Gb^—v) are functions that provide a cutoff a.t y = C(l), 
and ensure that a{x,t) and b{x,t) satisfy @ away from the reaction zone. The actual form 
of Ga^v) and GB^—y) is unimportant, provided that all moments of the tail for y > are 
defined and that there is no power-law behaviour for ?/ 0. 

This form leads to multiscaling behaviour for the moments of the particle profiles, by 
virtue of the power law tails of a and/or h when {m,n) 7^ (1, 1). Consider a function F of 
the form F(x,t) = (l){x/t°')G{x /t^^"^) , where 0(|/) as y — cx), — > 1 as y — > 0, 

G{y) ^ 1 as ?/ — > 0, and all positive moments of G are defined. Then the g'th moment of F 
is of the form 



1^°^ x^Fix, t) dx = xH~^(t) G (^-^ dx 



■ ^a{q+i)-s for/x>g + l, 

^a(g+i)-5iQg^ for/i = g + l. 



When yU > g + 1, the g'th moment of is finite, whereas, for /x < q + l, the dominant 
contribution comes from (x/t^^"^) = 0{1). Defining X^"^^ = [J x'^Fdx/ J F dxY^'^, we find 
that, for < 1, X^"^^ ~ t^/^ (with logarithmic corrections for /i = 1). For /i < 1, we find the 
multiscaling behaviour X^'^^ ~ f^i, with increasing monotonically as a function of q from 
a to 1/2. 

By substituting from ( p!^ the appropriate power-law tails of Ass and Bss, the multiscaling 
forms predict the following behaviour for the following quantities (without loss of generality, 
we have assumed m>n): 
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w 



J^oo x'^R{x, t)dx 
J^^R{x,t)dx 

2a 



t 

log t 

jy^a{x,t)dx 

\ogt 

(i-a)(A-l) 



for m — n < 3, 
for m — n = 3, 
for m — n > 3, 



(20) 



for 3m 
for 3m 



n < 5, 



n 



t 

\ogt 

t 



J^^x%{x,t)dx 
^^^h{x,t)dx 

t^" logt 

^l-(i-a)(K-l) 



for 3m >r2 + 5>m + 2, 
for m — n = 3, 
for m — n > 3, 



for m + 5 > 3n, 
for m + 5 = 3n, 
for m + 5 < 3n, 



(21) 



(22) 



/.2a 



SO 



where A = (ra + 2)/(m — 1), u = 2 + X, and k = {m + 2)/{n — 1). If these quantities were 
described by a one-length scaling theory, all of these quantities would behave as 
we describe departure from this behaviour as 'anomalous'. 
Defining 

JZoXR{x,t)dx 



coo 
-co ' 

fOO 
-OO ' 



JZ^R{x,t)dx' ^^^^ 

a similar procedure may be used to find the scaling behaviour for Xc- Notice, however, 
that the contribution to Xc coming from the scaling term is identically zero, since mRgg = 
Dd^a implies xRgg dx = 0. The behaviour of Xc is therefore determined wholly by the 
corrections to scaling. From (|l],^, one has 



xRdx = / — xdladx+ / — xd^bdx 
-OO J-oo m Jo n 



xdta dx 



xdfbdx. 



(24) 



The first two terms may be integrated by parts, yielding D[a{0,t)/m — b{0,t)/n\, which 
is zero by virtue of (|^) for the initial conditions u{0,0) = 0. The final two terms are of 
opposite sign, and cancel identically for m = n by symmetry. For m 7^ n, they typically 
have different scaling behaviours, and so the scaling behaviour is determined by the largest 
term. Substituting the forms 
scaling behaviour for Xc- 



T8|JT9|) for a and b, and differentiating, one finds the following 



Xr. 




for 2m — n < 4, 
for 2m — n = 4, 
for 2m — n > 4. 



(25) 
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III. NUMERICAL SIMULATIONS 



In order to verify the scaling ansatz (0-^ we have solved the react ion- diffusion equations 
(H,^ numerically. We approximated (|l]) and with a finite-difference method, with At, 
Ax satisfying kAt = 0.01, DAt/{Axf = 0.04, and a lattice of = 12001 sites, which was 
sufficient for finite-size effects to be unimportant. Henceforth we will choose the reference 
system in which Ax = At = 1, so that our numerical results will correspond to the case 
k = 0.01 and D = 0.04. The initial condition satisfied ao/m = bQ/n, so that Xf = from 
The reaction rate R thus obtained for times t = 10^ . . . 10'' and n = l, 2<m<4 are 
presented in Figures 1-3, rescaled by t""'"2 and plotted as a function of x/t"'. The solid curve 
in each figure is a numerical solution the steady-state equation ([l^), for the same values of 
D and k, and using the time dependent current as boundary condition. We see that the 
data appear to converge to the solid line as t — oo, so that the scaling ansatz d^-^) is valid 
in the sense that 

r+^R{x/r) F{x/e) as t ^ oo, (26) 

where F^x/f^) is a function of x/t". Notice that the convergence appears to be slower for 
larger values of m. 

Another interesting feature of Figures 1-3 is that the point at which R reaches its max- 
imal value differs from 0. As the x-axes of these figures have been rescaled by a factor t", 
the location of this point changes with time like t". This shows that the definition of x/ as 
the point of maximal reaction is not equivalent to the definition in this article as the point 
where u = 0. 

To verify the multiscaling properties of the reaction zone, we investigated the behaviour 
of the functions w, Wa, Wb, and Xc- According to (pO|-p^), they should diverge as t'^log'^it), 



with u and n being some exponents dependent on m and n. In Figs. 4-7 we have plotted 
these quantities, rescaled by log~'^t, on a log-log scale, using the theoretical values of k. 
The values of u^, Ua, I'b, and Uc (corresponding to the behaviour of the properties w, Wa, 
Wb, and Xc respectively) were estimated from a least-squares fit to the last decade in these 
figures, and are compared with the theoretical values in Table 1. Agreement for most of the 
exponents is very good. Since n = 1, the theory predicts that Ub assumes its non-anomalous 
value a, and the measured values are in close agreement. The values for z/^ are anomalous, 
and are also well reproduced. The values for i/^ and Uc are close to the theoretical predictions 
for m = 2, but deviate for m = 3 and 4. 

The reason why some of the exponents deviate from the theoretical values may be under- 
stood from Figs. 1-3. The convergence to the solid curve for x < 0, where Eq. (|T^ predicts 
an exponential decay, is very rapid — semilog plots of R{—x) were found to be in very good 
agreement with these predictions for up to 10 decades. However, the convergence of the pro- 
files to the steady-state profiles is much slower for x > (where the asymptotic behaviour 
is algebraic). In fact, numerical investigation of the steady-state profiles showed that the 
regime for which the power law behavior predicted by ([1^) appears is beyond the point at 
which the diffusive cutoff is already active in the data. This means that the asymptotic 
regime has not yet been reached. In view of this fact, the agreement between the measured 
exponents and the theory is surprisingly good. 



7 



IV. CONCLUSIONS 



The multiscaling theory predicts that the reaction profile of the system is described 
asymptotically by a scaling form, in the sense that Eq. (pB]) holds, but that the moments of 
the reaction and density profiles may have anomalous behavior. Numerical solutions of the 
react ion- diffusion equations verify the asymptotic scaling behavior, and also give values for 
the anomalous exponents close to those predicted by the theory. Longer times would have 
to be simulated to find better values for the exponents. 

The convergence to the asymptotic behavior becomes progressively slower as the order 
of the reaction is increased. This means that simulations probing the asymptotic behavior 
also become more difficult. Nevertheless, Eq. (^O]) suggests that it would be worthwhile to 
look at at least one case where m — n > 3. 



In the steady-state problem, it has been shown |Tl| that the critical dimension is dc = 
2/{m+n — l). This means that the reaction-diffusion equations correctly describe the scaling 
behavior of 'real' stochastic realizations for all physical dimensions, except for (m, n) = (2, 1) 
in dimension 1, where logarithmic corrections to the steady-state behavior are expected. 
Because of the strong link between the steady-state problem and the time-dependent problem 
studied in the present article, we expect the critical dimensions to be the same. The only 
microscopic simulation results available |^ agree broadly with the picture in the present 
article, though they are not of sufficiently high quality to give a thorough test. 

The theory outlined in this paper is based on heuristic arguments and numerical results 
only, so a first-principles analytical justification is needed. It would be possible to write 
the equation of motion for the corrections to the multiscaling terms, and then investigate 
whether they are truly small and do not contribute to the behavior of the moments. Although 
no conclusive results have been found, a preliminary attempt at this procedure suggests 
that there may be values for m and n where the assumption may break down that the 
interpenetration of the species is smaller than the number of particles that have reacted 



p!3| . The multiscaling theory presented in this paper might therefore have to be revised in 



those cases. 
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TABLES 



TABLE I. The measured values of the exponents i/^, Va and Vf, describing w, Wa, and Xc 
respectively, together with the non-anomalous exponent a, for three values of the duple {m,n). 
The values in square brackets are the predictions of the multiscaling theory; the presence of a value 
for K indicates the presence of a logarithmic correction of the form log''(i), account of which has 
been taken in calculating the i'. 



Quantity (2^1) (3^) (4^1) 

"4 10 3 

0.26 [i] 0.35 [^] 0.37 [i,K = i] 

Ua 0.27 [z,K = i] 0.44 [0.45] 0.50 [^,k--^^ 



i^b 0.25 [|] ' 0.29 [^] 0.32 [|i 

i^c -0.25 [-i] 0.06 [0] 0.21 [i] 
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FIGURES 

FIG. 1. R-t^ vs. x/t" forn = 1 and m = 2 and t = 10^, 10^ . . . , lO'^. 

FIG. 2. R-t'^ vs. x/i" forn = 1 and m = 3 and t = 10^, 10^, . . . , 10^. 

FIG. 3. R-tl^ vs. x/f^ forn = 1 and m = 4 and t = 10^, 10^, . . . , lO'^. 

FIG. 4. Log-log plot of wlog~'^{t) versus t. The straight lines are a least-squares fit to the 
last decade. 

FIG. 5. Log-log plot of ■Walog^'^{t) versus t. The straight lines are a least-squares fit to the 
last decade. 

FIG. 6. Log-log plot of wi,\og~'^{t) versus t. The straight lines are a least-squares fit to the 
last decade. 

FIG. 7. Log-log plot of Xclog~'*(t) versus t. The straight lines are a least-squares fit to the last 
decade. 
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